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Abstract 

CO ' 

This paper studies the newly emerging wireless powered communication network in which one hybrid access 

point (H-AP) with constant power supply coordinates the wireless energy/information transmissions to/from a set 

of distributed users that do not have other energy sources. A "harvest-then-transmit" protocol is proposed where all 

Q-T 

users first harvest the wireless energy broadcast by the H-AP in the downlink (DL) and then send their independent 

information to the H-AP in the uplink (UL) by time-division-multiple-access (TDMA). First, we study the sum- 
throughput maximization of all users by jointly optimizing the time allocation for the DL wireless power transfer 
versus the users' UL information transmissions given a total time constraint based on the users' DL and UL 
channels as well as their average harvested energy values. By applying convex optimization techniques, we obtain 
the closed-form expressions for the optimal time allocation to maximize the sum-throughput. Our solution reveals 
an interesting "doubly near-far" phenomenon due to both the DL and UL distance-dependent signal attenuation, 
where a far user from the H-AP, which receives less wireless energy than a nearer user in the DL, has to transmit 
J> . with more power in the UL for reliable communication. As a result, the maximum sum-throughput is achieved by 

allocating substantially more time to the near user than the far user, thus resulting in unfair rate allocation among 
users. To overcome this problem, we furthermore propose a new performance metric so-called common-throughput 

with the additional constraint that all users should be allocated with an equal rate regardless of their distances to 
"xl" . 

the H-AP. We present an efficient algorithm to solve the common-throughput maximization problem. Simulation 
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results demonstrate the effectiveness of the common-throughput approach for solving the uniquely new doubly 
near-far problem in wireless powered communication networks. 

Index Terms 

Wireless power, energy harvesting, throughput maximization, doubly near-far problem, TDMA, convex opti- 
mization. 

I. Introduction 

Traditionally, energy-constrained wireless networks, such as sensor networks, are powered by fixed 
energy sources, e.g. batteries, which have limited operation time. Although the lifetime of the network 
can be extended by replacing or recharging the batteries, it may be inconvenient, costly, dangerous (e.g., 
in a toxic environment) or even impossible (e.g., for sensors implanted in human bodies). As an alternative 
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solution to prolong the network's lifetime, energy harvesting has recently drawn significant interests since 
it potentially provides unlimited power supplies to wireless networks by scavenging energy from the 
environment. 

In particular, radio signals radiated by ambient transmitters is a viable new source for wireless energy 
harvesting. It has been reported that 3.5mW and luW of wireless power can be harvested from radio- 
frequency (RF) signals at distances of 0.6 and 11 meters, respectively, using Powercast RF energy-harvester 
operating at 915MHz [1J. Furthermore, recent advance in designing highly efficient rectifying antennas will 
enable more efficient wireless energy harvesting from RF signals in the near future [2J. It is worth noting 
that there has been recently a growing interest in studying wireless powered communication networks 
(WPCNs), where energy harvested from ambient RF signals is used to power wireless terminals in the 
network, e.g., [|3]-[|5]|. In 0, a wireless powered sensor network was investigated, where a mobile charging 
vehicle moving in the network is employed as the energy transmitter to wirelessly power the sensor 
nodes. In [|4]|, the wireless powered cellular network was studied in which dedicated power-beacons are 
deployed in the cellular network to charge mobile terminals. Moreover, the wireless powered cognitive 
radio network has been considered in [5|, where active primary users are utilized as energy transmitters 
for charging their nearby secondary users that are not allowed to transmit over the same channel due to 
strong interference. Furthermore, since radio signals carry energy as well as information at the same time, 
a joint investigation of simultaneous wireless information and power transfer (SWIPT) has recently drawn 
a significant attention (see e.g. [6J-[1 1 J and references therein). 

In this paper, we study a new type of WPCN as shown in Fig. [fl in which one hybrid access point (H- 
AP) with constant power supply (e.g. battery) coordinates the wireless energy/information transmissions 
to/from a set of distributed users that have no other energy sources. The channels between the H-AP and 
the users are assumed to follow the narrow-band block-fading model. All users are each equipped with 
a rechargeable battery and thus can harvest and store the wireless energy broadcast by the H-AP. Unlike 
prior works on SWIPT (|6l-[fTT|. which focused on the simultaneous energy and information transmissions 
to users in the downlink (DL), in this paper we consider a different setup where the H-AP broadcasts 
only wireless energy to users in the DL while the users transmit their independent information using their 
individually harvested energy to the H-AP in the uplink (UL). We are interested in maximizing the UL 
throughput of the aforementioned WPCN by optimally allocating the time for the DL wireless energy 
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Fig. 1. A wireless powered communication network (WPCN) with wireless energy transfer (WET) in the downlink (DL) and wireless 
information transmissions (WITs) in the uplink (UL). 



transfer (WET) by the H-AP and the UL wireless information transmission (WIT) by the users. 
The main contributions of this paper are summarized as follows: 

• We propose a protocol termed "harvest-then-transmit" for the WPCN as shown in Fig. [Q where 
the H-AP first broadcasts wireless energy to users in the DL, and then the users transmit their 
independent information to the H-AP in the UL using their individually harvested energy by time- 
division-multiple-access (TDMA). 

• With the proposed protocol, we first maximize the sum-throughput of the WPCN by jointly optimizing 
the time allocated to the DL WET and the UL WITs given a total time constraint, based on the users' 
DL and UL channels as well as their average harvested energy amount. It is shown that the sum- 
throughput maximization problem is convex, and therefore we derive closed-form expressions for the 
optimal time allocation by applying convex optimization techniques Ifl2ll . 

• Our solution reveals a uniquely new "doubly near-far" phenomenon in the WPCN, which occurs 
when a far user from the H-AP which receives less amount of wireless energy than a near user in 
the DL has to transmit with more power in the UL for achieving the same information rate due to 
the doubly distance-dependent signal attenuation in both the DL WET and UL WIT. Consequently, 
the sum-throughput maximization solution allocates substantially larger time to the near user than 



the far user, thus resulting in unfair achievable rates among users. 
• To overcome the doubly near- far problem, we furthermore propose a new performance metric referred 
to as common-throughput with the additional constraint that all users should be allocated with an 
equal rate in their UL WITs regardless of their distances to the H-AP. We propose an efficient 
algorithm to maximize the common-throughput of the WPCN by re-optimizing the time allocated 
for the DL WET and UL WITs. By comparing the maximum sum-throughput versus the maximum 
common-throughput, we characterize the fundamental throughput-fairness trade-offs in the WPCN. 
The rest of this paper is organized as follows. Section [I] presents the WPCN model and the pro- 
posed harvest-then-transmit protocol. Section [III] studies the sum-throughput maximization problem, and 
introduces the doubly near- far phenomenon. Section HVl formulates the common-throughput maximization 
problem and presents an efficient algorithm to solve it. Section |V] presents simulation results on the 
sum-throughput versus common-throughput comparison. Finally, Section |VT] concludes the paper. 

II. System Model 

As shown in Fig. [fl this paper considers a WPCN with WET in the DL and WITs in the UL. The 
network consists of one H-AP and K users (e.g., sensors) denoted by Uu % = 1, • • • , K. It is assumed 
that the H-AP and all user terminals are equipped with one single antenna each. It is further assumed 
that the H-AP and all the users operate over the same frequency band. In addition, all user terminals are 
assumed to have no embedded energy sources; thus, the users need to replenish energy from the received 
signals broadcast by the H-AP in the DL, which is then used to power operating circuits and transmit 
information in the UL. 

The DL channel from the H-AP to user C/j and the corresponding reversed UL channel are denoted 
by complex random variables hi and g it respectively, with channel power gains hi = \hi\ 2 and g^ = \cn\ 2 . 
It is assumed that both the DL and UL channels are quasi-static flat-fading, where hi's and gi's remain 
constant during each block transmission time, denoted by T, but can vary from one block to another. It 
is further assumed that the H-AP knows both hi and g„ i — 1, ■ ■ • , K, perfectly in each block. 

The network adopts a harvest-then-transmit protocol as shown in Fig. |2] In each block, the first r T 
amount of time, < t < 1, is assigned to the DL for the H-AP to broadcast wireless energy to all users, 
while the remaining time in the same block is assigned to the UL for information transmissions, in which 
users transmit their independent information to the H-AP by TDMA. The amount of time assigned to user 
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Fig. 2. The harvest-then-transmit protocol. 



Ui in the UL is denoted by r(T, < Tj < 1, i — 1, • • - K. Since r , T\, • ■ ■ , Tk represent the time portions 
in each block allocated to the H-AP and users U\, ■ ■ ■ , Uk for UL WET and DL WITs, respectively, we 
have 

K 

J><1. (1) 

i=0 

For convenience, we assume T = 1 in the sequel without loss of generality. 

During the DL phase, the transmitted complex baseband signal of the H-AP in one block of interest 
is denoted by xa- We assume that xa is an arbitrary random signal^ satisfying E[|xa| ] = Pa, where Pa 
denotes the average transmit power at the H-AP The received signal at Ui is then expressed as 



Vi 



hiXA T Zi 



, K, 



(2) 



where yi and z% denote the received signal and noise at Ui, respectively. It is assumed that Pa is sufficiently 
large such that the energy harvested due to the receiver noise is negligible. Thus, the amount of energy 
harvested by each user in the DL can be expressed as (assuming unit block time, i.e., T = 1) 



Ei = (iP A hiT , i = 1, 



K. 



(3) 



where < Q < 1, i = 1 • • • , K, is the energy harvesting efficiency at each user. For convenience, it is 
assumed that d = ■ ■ ■ = (k = C m me sequel of this paper. 

After the users replenish their energy during the DL phase, in the subsequent UL phase they transmit 
their independent information to the H-AP in their allocated time slots. It is assumed that at each user 
terminal, a fixed portion of the harvested energy given by © is used for its information transmission in the 
UL, denoted by r\i for Ui, < rji < 1, i = 1, ■ ■ • , K. Within r* amount of time assigned to Ui, we denote 
Xi as the complex baseband signal transmitted by U if i — 1, ■ • • , K. We assume Gaussian inputs, i.e., 

'Note that x a can also be used to send DL information at the same time; however, this usage will not be considered in this paper. 
Interested readers may refer to recent works on SWIPT llol- lll II . 



Xi ~ CAf(0,Pi), where CJ\f(fj,,a 2 ) stands for a circularly symmetric complex Gaussian (CSCG) random 
variable with mean /x and variance a 2 , and Pi denotes the average transmit power at Ui, which is given 
by 

71; Ea 

Pi=—, i = l, ■■■, K. (4) 

n 

For the purpose of exposition, we assume r\i = 1, Vi, in the sequel, i.e., all the energy harvested at each 
user is used for its UL information transmission. The received signal at the H-AP in the zth UL slot is 
then expressed as 

UA,i = y/OiXi + z A , % = 1, • • • , K, (5) 

where y Aji and z A denote the received signal and noise at the H-AP, respectively. It is assumed that z A ~ 
CA/"(0, a 2 ). From ©-©, the achievable throughput of Ui in bits/second/Hz (bps/Hz) can be expressed as 

9i Pi 



Ri (t) = nlogo 1 + 
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Ta 2 



= T;log 2 ( l + 7<— j , i = l,---,K, (6) 

where r = [t T\ ■ ■ -r K ], and T represents the signal-to-noise ratio (SNR) gap from the additive white 
Gaussian noise (AWGN) channel capacity due to a practical modulation and coding scheme (MCS) used. 

In addition, % is given by 

Chi9iP A . v . 

li = Ta2 , i = l, ••• , K. (7) 

From ©, it is observed that Ri (r) increases with r for a given r». In addition, it can also be shown 
that Ri (t) increases with r, for a given t . However, r and r»'s cannot increase at the same time given 
their sum-constraint in (OQ). Fig. |3] shows the throughput given in © for the special case of one single 
user in the network, i.e., K = 1, versus the time allocated to the DL, r , with 71 = lOdB, assuming that 
(OQ) holds with equality. It is observed that the throughput is zero when r = 0, i.e., no time is assigned 
for WET to the user in the DL and thus no energy is available for WIT in the UL, as well as when 
r = 1 or n = 1 — t = 0, i.e., no time is assigned to the user for WIT in the UL. It is also observed 
that the throughput first increases with To when r < Tq = 0.42, but decreases with increasing t when 
r > Tq, where Tq is the optimal time allocation to maximize the throughput. This can be explained as 
follows. With small r , the amount of energy harvested by U\ in the DL is small. In this regime, as U\ 
harvests more energy with increasing r , i.e., more energy is available for the information transmission 
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Fig. 3. Throughput versus time allocated to DL WET in a single-user WPCN with 71 = lOdB. 



in the UL, the throughput increases with t . However, as r becomes larger than Tq, the throughput is 
decreased more significantly due to the reduction in the allocated UL transmission time, T\ = 1 — r ; as 
a result, the throughput starts to decrease with increasing r . Therefore, there exists a unique optimal Tq 
to maximize the throughput. 

III. Sum-Throughput Maximization 
In this section, we maximize the sum-throughput of the WPCN modeled in Section HH From ©, the 

K 

sum-throughput of all users is given by i? sum (t) = ^2 Ri (r), which is a function of the DL and UL 

i=l 

time allocation r. Therefore, from (H) the sum-throughput maximization problem is formulated as 



(PI) 



max R sum (t) 

T 

K 

s.t. 7~Vj < 1, 

i=0 



r, > 0, 2 = 0, 1, 



• if. 



(8) 



Lemma 3.1: Ri (r) is a concave function of r for any given i G {1, 
Proof: Please refer to Appendix [A] 



*}• 




5 4 4.5 5 



Fig. 4. Plot of / (z) given in l[9} versus z > 0. 



From Lemma 13. 1 L it follows that R sum (r ) is also a concave function of r since it is the summation of 
Ri (t)'s. Therefore, (PI) is a convex optimization problem, and thus can be solved by convex optimization 
techniques. To solve (PI), we first have the following lemma. 

Lemma 3.2: Given A > 0, there exists a unique z* > 1 that is the solution of / (z) = A, where 



f{z) = z\nz-z + l, z>0. 



(9) 



Proof: Please refer to Appendix |B] ■ 

Fig. H] shows / (z) given in © with z > 0. It is observed that / (z) is a convex function over z > 

where the minimum is attained at z = 1 with / (1) = 0. Therefore, given < A < 1, there are two 

solutions for / (2) = A, among which one is smaller than 1 and the other is larger than 1, i.e., z* > 1. 

On the other hand, if A > 1, there is only one solution for / (z) = A, which is larger than 1, i.e., z* > 1. 

The above observations are thus in accordance with Lemma 13.21 

Proposition 3.1: The optimal time allocation solution for (PI), denoted by r* = [tq t* • 
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given by 
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i = Q 



(10) 



K 



where A = Y^ li > and z* > 1 is the corresponding solution of f (z) = A as given by Lemma 13.21 



i=l 



Proof: Please refer to Appendix | 



It is worth noting that A > always holds since from © we have j t > 0, i — 1, 2, • • • , K, provided 
that hi 7^ and gi ^ 0. Hence, given a set of strictly positive 7/s, according to Lemma 13.21 z* > 1 is 
uniquely determined with the presumed A, thus resulting in a unique solution t* for (PI), with r* > 0, 
2 = 0, 1, • • • , K, i.e., the time allocated to the DL WET is always greater than zero, and so is the time 
allocated to each user in the UL WIT. Furthermore, from Proposition [3J] we have the following corollary. 

Corollary 3.1: In the optimal time allocation solution of (PI), Tq is a monotonically decreasing function 
of A > 0. 

Proof: Please refer to Appendix [D] ■ 

From Corollary 13.11 it is inferred that the time allocated to the DL WET decreases with increasing 

K 

7i's, or channel power gains h^s and/or g/s, since A = ^ 7^ and ji oc hiQi, % — 1, • • • , K, as shown in 
©. As a result, r*'s, z = 1, ■ ■ • , -ft', increase with A, i.e., the time allocated to the UL WIT increases 
with 7j's. This is an interesting observation implying that when the channel power gains, /ij's and g^s, 
become larger, we should allocate more time to the UL WITs instead of the DL WET to maximize the 
sum-throughput. This is because with larger Tj's, the required energy for UL WITs becomes smaller given 
any transmission rate; thus, each user can harvest sufficient amount of wireless energy from the H-AP 
even with a smaller time allocated to the DL WET. 

Note that the sum-throughput maximization solution given in (flOl) allocates more time to a 'near' user 
to the H-AP than a 'far' user, since hi oc D^ ad , g { oc D^ au , and 7, oc /i^ in ©, where a^ > 2 and 
a u > 2 denote the channel pathloss exponents in the DL and UL, respectively, and D ri denotes the distance 
between the H-AP and U l . Thus, from CEO]) it follows that t* oc D^ ad+au) , i = 1, ■ ■ ■ , K, which results 
in an unfair time and rate allocation among users in a WPCN, a phenomenon termed "doubly near-far 
problem". To further illustrate this problem, Fig. |5] shows the sum-throughput R sum (r) versus UL WIT 
time allocation n and r 2 for a two-user network with K = 2 and Di = \D 2 . It is assumed that the channel 
reciprocity holds for the DL and UL and thus hi = gi, i = 1, 2, with ad = a u = 2. Accordingly, we set 
71 = 22dB and 72 = lOdB, with 71/72 = {D l /D 2 ) ad+au . It is observed from Fig. |5] that R snm (r) = 
when 7"i = r 2 = 0, t = 1 — (r x + r 2 ) = 1, or T\ + r 2 = 1, r = 0, since no time is allocated to the users 
for the UL WITs in the former case, while no time is allocated to the DL WET in the latter case. The 
numerical result of sum-throughput clearly shows that -R S um (t) is strictly positive when < T\ + r 2 < 1. 
In addition, it is observed that the optimal DL and UL time allocation to maximize the sum-throughput 
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Fig. 5. Sum-throughput versus time allocation. 



is r* = [0.2441, 0.7114, 0.0445] where rf = 16r 2 \ i.e., r* = (D 2 /D 1 ) a * +a <"r$, which is consistent 
with (flOT) . Furthermore, at the optimal r*, i? sum (r*) = 4.58 bps/Hz, with i? x (r*) = 4.13 bps/Hz and 
-R2 (t*) = 0.45 bps/Hz, which demonstrates the very unfair rate allocation between the two users due to 
the doubly near-far problem. 

For comparison, consider UL transmissions in a conventional TDMA-based wireless network for WIT 
only lfT3l - lfi~6ll . where each user is equipped with a constant energy supply, and has a given energy 
consumption at each block denoted by E. It then follows from © that 7$ oc D^ au and thus from (flOT ). the 
optimal time allocation to maximize the sum-throughput of such a conventional TDMA network should 
satisfy t* oc D~ au , i = 1, • • • , K, and r * = (since no DL WET is needed). Clearly, the WPCN suffers 
from a more severe near- far problem than the conventional TDMA network. With the same setup as for 
Fig. [51 in Fig. |6] we show the optimal throughput of U 2 , R2 ( r *)> normalized by that of U\, R\ (r*), in 
a WPCN versus that in a conventional TDMA network for different values of the pathloss exponent a, 
with ad = a u = a. For the WPCN, 71 is set to be fixed as 71 = 22dB for Ui, while 72 = 10, 7, 4, 1, 
and — 2dB for a = 2, 2.5, 3, 3.5, and 4, respectively, since 71/72 = (D 2 /Di) a . For the conventional 
TDMA network, E for both U\ and U 2 are assumed to be E = ^(E 1 + E 2 ) with E 1 and E 2 denoting the 
average harvested energy at XJ\ and U 2 in the WPCN under comparison, respectively; it then follows that 
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Fig. 6. Throughput ratio in a two-user WPCN versus conventional TDMA network with constant energy supply. 



for the conventional TDMA network, 71 = 13dB and 7 2 = 7.1, 5.5, 4.0, 2.4, and 0.9dB for a = 2, 2.5, 
3, 3.5, and 4, respectively. From Fig. |6] it is observed that the rate ratio of the two users in the WPNC 
case decreases twice faster than that in the conventional TDMA in the logarithm scale due to the more 
sever (doubly) near-far problem. 

IV. Common-Throughput Maximization 

In this section, we tackle the doubly near-far problem in the WPNC by applying the common-throughput 
maximization approach, which guarantees an equal rate allocation to all users. From (Q3 and ©, the 
common-throughput maximization problem is formulated as 



(P2) : max R 



s.t. Ri(T)>R, i = l---K, 



(ID 



r ev, 

where R denotes the common-throughput and V is the feasible set of r specified by (OQ) and ([8]). It can 
be easily shown that the optimal time allocation solution r* for (P2) should allocate the same optimal 

K 

throughput to all the users, denoted by R* = R 1 (r*) = • ■ • = Rk (f *), with J2 r* = 1. 

i=0 



12 

The maximum common-throughput R* is the maximum of all the feasible common-throughput R that 
satisfies the rate inequalities in (fTTT) of (P2). To solve (P2), given any R > 0, we first consider the 
following feasibility problem: 

Find r 

s.t. Ri(r)>R, i = l, ■•■ ,K, 

TEV. (12) 

Since the problem in (fT2l is convex, we consider its Lagrangian given by 

K 

C(t,\) = -Y,\(Ri{T) ~ R), (13) 

where A = [Ai, • • • , \k\ > ('here >' denotes the component- wise inequality) consists of the Lagrange 



multipliers associated with the K rate constraints in (|T2l . The dual function of problem ([121 1 is then given 

by 

£(A) = min£(r,A). (14) 

The dual function Q (A) can be used to determine whether problem (fT2l is feasible, as provided in the 
following lemma. 

Lemma 4.1: For a given R > 0, problem (fT2l) is infeasible if and only if there exists an A > such 
that g(X) > 0. 

Proof: Please refer to Appendix |E] ■ 

Next, we obtain G(X) in (fl4l) for a given A > by solving the following weighted sum-throughput 
maximization problem, which follows from ([TJt . 



K 

max 

T 



^A^(t) 

S.t. T G P. (15) 



Like (PI), the problem in (1151) is convex and thus can be solved by convex optimization techniques. 
Similar to Proposition 13.11 for the sum-throughput maximization case with \ = 1, Vz, we obtain the 
optimal time allocation solution for the weighted sum-throughput maximization problem in (fl3T >. given in 
the following proposition. 
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Proposition 4.1: Given A > 0, the optimal time allocation solution for (fT5T ), denoted by t* = [tq, t*, 



1 + E (7i/4) 

3=1 



li/ Z t 






if, (17) 



where zf,i — l, ■ ■ ■ , If, is the solution of the following equations: 



In(l + ^)--5_ = ^-ln2 ) (18) 

1 + Zj Aj 



A' 



Aj7-. 



E^t = "*^. a» 



with fi* > being a constant. 

Proof: Please refer to Appendix IB ■ 

With Proposition 14. 1[ we can compute t* efficiently as shown next. Denote the left-hand sides (LHSs) 
of (fT8~l) and (fT9b as Q (zi) and S (z) with z = [24, z 2 , ••• , £«"]> respectively. Note that Q (zi) is an 
increasing function of Z{, % — 1, • • ■ , if (see Appendix 0, whereas 5 (z) is a decreasing function with 
respect to each individual Z{. Given any fi > 0, suppose that Zi is the solution of Q (zi) = f-hi2, 
i = 1, ■ ■ • , K , in (fTBl . With these Zj's, there are two possible cases to consider as follows. If in (fl9l) 
the resulting S (z) > /iln2, we should increase /i since Zj's satisfying Q (zi) = #-ln2, i = 1, • • • , if , 
will increase with // given that Q (<%), Vi, is an increasing function of zf, as a result, S 1 (z) will decrease 
since it is a decreasing function of each individual Z{. Otherwise, /i should be decreased to satisfy (fT9l ) if 
S (z) < /iln2. Therefore, z*'s and //* can be obtained by iteratively updating z/s and /i as above until 
convergence is reached. Then, r* can be computed from (fTol and (TTvT > accordingly. 

Given R, A, and the obtained r* by solving the problem (TT3T ) with Proposition 14.11 we can compute 
the corresponding i?j (r*), i = 1, • • • , if, and thus £ (A) in (fT4l) using (TT3T) . If C? (A) > 0, it follows 
from Lemma |4~T1 that problem (fT2l is infeasible, i.e., R > R*. Therefore, we should decrease R and 
solve the feasibility problem in (fT2l) again. On the other hand, if C? (A) < 0, we can update A using 
sub-gradient based algorithms, e.g. the ellipsoid method [fTTll . with the subgradient of Q (A), denoted by 
v = \v\ v-i • • ■ vk] t , with 

Vi = r*log 2 (l + ^pj -R, 1 < i < if, (20) 



I'l 



TABLE I 

Algorithm to solve (P2). 



r 



1) Initialize i? min = 0, _R max > R* 

2) Repeat 

2. Initialize A > 0. 

3. Given A, solve the problem in (Q3} by Proposition 14. 1 1 

4. Compute Q (A) using (13). 

5. If 5 (A) > 0, R is infeasible, set i? ma x «— -R, go to step 1. 
Otherwise, update A using the ellipsoid method and the 
subgradient of Q (A) given by ( |20t . If the stopping criteria 
of the ellipsoid method is not met, go to step 3. 

6. Set R mln <r- R. 



3) Until R n 



■ Rmln < S, where 5 > is a given error tolerance. 



until A converges to A* with A* denoting the maximizer of Q (A) or the optimal dual solution for the 
problem (fl"2l . If Q (A*) < 0, it then follows that the problem in (fl"2l is feasible and R < R*. In this case, 
R should be increased for solving the feasibility problem in (TT2|) again. Consequently, R* can be obtained 
numerically by iteratively updating R by a simple bisection search |fl2|. To summarize, one algorithm to 
solve (P2) is given in Table HI 

Fig. [7] shows the common-throughput in bps/Hz versus n and r 2 for the same two-user channel setup as 
for Fig. |U It is observed that the optimal time allocation for (P2) is given by r* = [0.3683, 0.1386, 0.4932], 
which results in R* = Ri (r*) = i?2 (t*) = 1.46bps/Hz. Comparing to Fig. |5] where the sum-throughput is 
maximized, the time portion allocated to the near user, U%, is decreased substantially from 0.7114 to 0.1737, 
while that to the far user, U 2 , is greatly increased from 0.0445 to 0.4669. Consequently, the throughput 
of U 2 is increased from 0.45bps/Hz to 1.46bps/Hz, while that of U\ is decreased from 4.13bps/Hz to 
1.46bps/Hz, the same rate as U\. This result shows the effectiveness of the proposed common-throughput 
approach to tackling the doubly near-far problem in a WPCN. 

Fig. [8] shows the ratio of the optimal time allocated to the far user U 2 over that to the near user U\, 
i.e., t 2 /t*, in (PI) versus (P2), with different values of the common pathloss exponent a in both DL 
and UL, where the same two-user channel setup as for Fig. |5]is considered. It is assumed that 71 for the 
near user U\ is fixed as ji = 22dB and 72 for the far user U 2 is set as for Fig. |6] It is observed that 
the time ratio of (PI) in the logarithm scale decreases linearly with a to maximize the sum-throughput, 

2 The initial value of -R max can be chosen as any arbitrary large number such that it satisfies 7? m ax > R*. 
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Fig. 7. Common-throughput versus time allocation. 
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Fig. 9. Throughput region of a two-user WPCN with (a) corresponding to the maximum sum-throughput; and (b) corresponding to the 
maximum common-throughput. 



which can also be inferred from (fTOl ). due to the doubly near-far problem. In contrary, t^/t* is observed 
to increase with a to maximize the common-throughput in (P2), since more time is allocated to the far 
user f/ 2 instead of the near user U\ as the ratio between 71 and 72, i.e., 71/72, increases with a. 

Notice that (PI) and (P2) deal with two extreme cases of rate allocation to the users where the fairness 
is completely ignored and a strict equal fairness is imposed, respectively. More generally, Fig. |9] shows the 
achievable throughput region of a two-user WPCN by solving the weighted sum-throughput maximization 
problem in (fT5l) with different rate weights for the near and far users, under the same channel setup as for 
Fig. [5] It is observed that the boundary of the throughput region characterizes all the optimal throughput- 
fairness trade-offs in this two-user WPCN, which include the rate pairs obtained by solving (PI) for the 
maximum sum-throughput as well as by solving (P2) for the maximum common-throughput, shown as 
points (a) and (b) in the figure, respectively. It is worth noting that the common-throughput approach 
for characterizing the achievable rate region of multi-user communication systems under strict fairness 
constraints can be considered as one special case of the "rate-profile" method as proposed in [18J. 
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Fig. 10. Sum-throughput vs. common-throughput. 



V. Simulation Result 



In this section, we compare the maximum sum-throughput by (PI) versus the maximum common- 
throughput by (P2) in an example WPCN with K = 2, D\ = 5m, and D 2 = 10m. The bandwidth is set 
as 1MHz. It is assumed that the channel reciprocity holds for the DL and UL and thus hi — gu i — 1, 2, 
with the same pathloss exponent a d = a u = a. Accordingly, both the DL and UL channel power gains 
are modeled as hi = gi = 1CT 2 p\D~[ a , % = 1,2, where p, ; represents the additional channel short-term 
fading which is assumed to be Rayleigh distributed, and thus p\ is an exponentially distributed random 
variable with unit mean. Note that in the above channel model, a 20dB average signal power attenuation 
is assumed at a reference distance of lm. The AWGN at the H-AP receiver is assumed to have a white 
power spectral density of — 160dBm/Hz. For each of two users, the energy harvesting efficiency for WET 
is assumed to be C = 0.5. Finally, we set T = 9.8dB assuming that an uncoded quadrature amplitude 
modulation (QAM) is employed. 

Fig. \W\ shows the maximum sum-throughput versus the maximum common-throughput in the same 
WPCN for different values of transmit power at H-AP, Pa, in dBm, by averaging over 1000 randomly 
generated fading channel realizations, with fixed a = 2. As shown in Fig. \W[ when the sum-throughput 
is maximized, the throughput of U\ dominates over that of U 2 due to the doubly near-far problem, which 
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Fig. 11. Throughput vs. pathloss exponent. 



results in notably unfair rate allocation between the near user (<7i) and far user (U 2 ) in this example. It 
is also observed that the maximum common-throughput for the two users is smaller than the normalized 
maximum sum-throughput (to the number of users, K = 2), which is a cost to pay in order to ensure a 
strictly fair rate allocation to the two users regardless of their distances from the H-AP. 

Next, by fixing Pa = 20dBm, Fig. \TT\ shows the throughput comparison for different values of the 
common pathloss exponent a in both the DL and UL. It is observed that when the sum-throughput is 
maximized, the throughput of the near user U\ converges the maximum sum-throughput as a increases, 
whereas that of the far user <7 2 converges to zero, which indicates that the WPCN suffers from a more 
severe unfair rate allocation between the near and far users as the pathloss exponent increases, due to the 
doubly near-far problem. In addition, the maximum common-throughput for the two users is observed to 
decrease faster with increasing a than the normalized maximum sum-throughput. This is because as a 
increases, (P2) allocates more time to the far user U 2 instead of near user V\ in order to ensure the equal 
rate allocation among users since the ratio 71/72 increases with a, whereas (PI) allocates more time to 
V\ instead of <7 2 as a increases. 
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VI. Conclusion 

This paper studied a new type of wireless RF (radio frequency) powered communication network with a 
harvest-then-transmit protocol, where the H-AP first broadcasts wireless energy to distributed users in the 
downlink and then the users transmit their independent information to the H-AP in the uplink by TDMA. 
We revealed an interesting new phenomenon in such hybrid networks, so-called doubly near-far problem, 
which occurs when a far user from the H-AP with less wireless energy harvested than a near user has to 
transmit with more power for reliable communication with the same rate due to the distance-dependent 
signal attenuation in both the downlink WET and uplink WIT. As a result, notably unfair time and 
throughput allocation among the users occurs when the conventional metric of network sum-throughput 
is maximized. To overcome this problem, we proposed a new common-throughput approach to allocate 
equal rates to all users regardless of their distances from the H-AP by allocating the transmission time 
to users inversely related to their distances to the H-AP. Simulation results showed that this approach is 
effective in solving the doubly near-far problem in the WPCN, but at a cost of certain sum-throughput 
degradation. 

Appendix A 
Proof of Lemma I3TT1 

Denote the Hessian of Ri (r) defined in © as 

0<j,m<K, 



V 2 R % (r) 



d {i) 



where d^ m denotes the element of V 2 Ri (r) at the jth row and mth column. From ©, the diagonal 
entries of V 2 i?j (r), i.e., j = m, can be expressed as 

4:1 = -I^7^ 3 /3r 2 , j = i (21) 

I , otherwiese, 

where /% = 1 + 2£21 . In addition, the off-diagonal entries of V 2 i?j (t) can be expressed as 

,(0 ,(i) J b^ToTf 2 /?r 2 , 3 = i and m = 

u 7,m — a m.j ~ ) _ , , V^J 

| (J , otherwise. 

Given an arbitrary real vector v = [v , v%, ■ ■ ■ , vk} t , since t, > 0, it can be shown from (12TT ) and 
(1221) that 
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v V ^ (r) V = -j^^2 ^3 (^ r o - ^OTi) 

<o, 

i.e., V 2 i?j (r), Vi, is a negative semidefinite matrix. Therefore, Ri (r) is a concave function of r 
[r T\ ■ ■ ■ r N ] T [1T21 . This completes the proof of Lemma [3TT1 

Appendix B 
Proof of Lemma [372] 



From (|9]), we have 



lim/(z) = l, (23) 

^ = m 2 , (24) 

^W- 1 (25) 



<9z 2 z 

Thus it follows from (|24l ) and (12~5T ) that / (z) is a convex function over z > and its minimum is attained 
at z = 1 with / (1) = 0. This implies that / (z) > with 2 > 1 and is monotonically increasing with z 
in this regime. Therefore, given A > 0, / (z) — A has a unique solution z* > 1. This completes the proof 
of Lemma 13.21 

Appendix C 
Proof of Proposition |3JJ 

The Lagrangian of (PI) is given by 

4™ (t, i/) = i?sum (r) -v I yVj- 1 J , (26) 

where i/ > denotes the Lagrange multiplier associated with the constraint in (Q]). The dual function of 
(PI) is thus given by 

g (u) = min £ sum (r, v) , (27) 

where V is the feasible set of r specified by ([T]) and ©. 

It can be shown from ([T]) and d8]) that there exists an t E T> with Tj > 0, z = 0, 1, • • • K, satisfying 

K 

Y^ Ti < 1, and thus strong duality holds for this problem thanks to the Slater's condition [fl2l . Since (PI) 

2=0 
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is a convex optimization problem for which the strong duality holds, the Karush-Kuhn-Tucker (KKT) 
conditions are both necessary and sufficient for the global optimality of (PI), which are given by 



A' 



X>;<i> (28) 

t=0 

i>;-l)=0, (29) 

8=0 / 

^ sum (r*) -u* = 0, % = 0, ■ • • , K, (30) 



t=0 

'(■' 

\t=0 

d_ 

Or, 



where r*'s and v* denote the optimal primal and dual solutions of (PI), respectively. It can be easily 

K 

verified that J^t* = 1 must hold for (PI) and thus from (|29l without loss of generality, we assume 

i=0 

v* > 0. From ©, it follows that 

K 

V 7 ' r . =t/*ln2, (31) 

tTl + 7^ 



M7i4 = ^* ln2 , 1<^<^, (32) 



where t (x) is defined as 



t(x) = \n(l+x) — , x>0. (33) 

1 + x 



Given 1 < i,j < K, from (1311 we have 



Try \ -73 



M7«i =M7i4 > ^i- ( 34 > 



It can be easily shown that £ (x) is a monotonically increasing function of x > since <i£ (a;) /da; = 
x(l + x)~ > for x > 0. Therefore, equality in (|34l ) holds if and only if 7^ = 7./-%, 1 < i,j < K i.e., 

*=£ = .. .2*=C. (35) 

' 1 r 2 'if 

if 

Note that 1 — Tq = J^ r* and r* = — r* from (|28l ) and (1331 ), respectively. Therefore, r* can be expressed 

as 

< = (1 - r *) -^- = (1 - r*) * (36) 

E7i 

K 

where A = J^ 7,. In addition, it follows from (|3Q|) . (|35T) . and (|36l) that 

i=i 

In (1 + Ct*) -£— = — - , Ki< K, (37) 

V 0J 1 + CV * 1 + CV*' ~ ~ 
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where C is defined in (T3~5b . Since C = -r-^-r from (T3~6b , we can modify (T3"7b as 

zlnz-z- A+l = 0, (38) 



where z = 1 + -r-V It is observed that z > 1 if A > and < Tn < 1. From Lemma [3T2l there exists a 

1— r„ u 

unique z* > 1 that is the solution of (1381) . Therefore, the optimal time allocation to the DL WET is given 

by 

z* — 1 

To = -A — ■ (39) 

A + z* - 1 

In addition, from (l36l) and (139b , the optimal time allocation to the UL WITs r*, 1 < i < K, is given by 

T - = ATT^l- <40) 

This thus proves Proposition 13.11 

Appendix D 
Proof of Corollary 13.11 

It can be easily shown from (flOb that Tq — 1 with A = 0. From (|38b and (|39b , Tq can be alternatively 
expressed as 

To = ^^- (41) 

2*mz* 



Given A > and thus z* > 1, both z* lnz* and z* — 1 in (|41b increase with A since z* increases with A 
as shown in the proof of Lemma [3T2l given in Appendix El Furthermore, since 4-(zlnz) = 1 + \nz and 
■£(z — 1) = 1, it follows that -£^(z*\nz*) > -J^{z* — 1) with z* > 1, i.e., z*\nz* increases faster with 
z* than z* — 1. Therefore, it can be verified that z* In z* increases faster with A than z* — 1, and thus Tq 
decreases monotonically with increasing A. Finally, it can be shown that Tq — > as A — > oo from the 
fact that z* increases with A. 

This thus completes the proof of Corollary 13.11 

Appendix E 
Proof of Lemma |4J] 

We first prove the "if" part of Lemma |4~T1 If r' G V is a feasible solution for ([Tib given R > 0, i.e., 
-Rj (t') > -R, i = 1, ■ ■ • , K, then for any A > it follows from CCD that 

£(A)<£(t',A)<0, 
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and thus max Q (A) < 0, which contradicts with the given assumption that there exists an A > such 

A>0 

that Q (A) > 0. The "if" part is thus proved. 

Next, we prove the "only if" part of Lemma 14.11 by showing that its transposition is true, i.e, the 
problem in ([T2T ) is feasible if Q (A) < 0, VA > 0, by contradiction. Suppose that problem (|T2T> is feasible 
and there exits an A" > where Q (A") > 0. However, since (PT21) is assumed to be feasible, there exists 
an t" e V such that R { (r") > R, Vz, resulting in A/' fa (r") - R) > since A" > 0. From (H3]) and 
(PT41) . we thus have 

Q (A") < - Y^U A *" (^ ( r ) - ^) < o. 

This contradicts Q (A") > 0, and thus problem (jH) is feasible if Q (A) < 0, VA > 0. The "only if" part 
is thus proved. 

Combining the above proofs of both "if and "only if parts, Lemma |4~T1 thus follows. 

Appendix F 
Proof of Proposition |4JJ 

Given A > 0, the Lagrangian of problem (fl"5T ) is given by 

K /K \ 

£ W SR ( r > Z 1 ) = Y X%Ri ( T ) ~ ^\Y Ti ~ 1 ) ' ^ AT> 

i=l \i=0 / 

where fi > denotes the Lagrange multiplier associated with the constraint in ©. The dual function of 



problem (1151) is thus given by 

£wsr O) = min £ W sr (t, //) . (43) 

Similar to (PI), it can be easily shown that the problem in (fT5l ) is convex with zero duality gap. Therefore, 
the following KKT conditions must be satisfied by the optimal primal and dual solutions of problem <TT3T >: 

In 1 + ft-J V = V ln2 ' i = h---,K, (44) 

V rt) 1 + Pi % A, 

f;-^r = Ai*ln2, (45) 

i=0 

where /i* > is the optimal dual solution. We then obtain (fT8l) and (TT91) by changing variables as 

Zi = H% i = l, ■■■ ,K, (47) 

r 



24 



in (|44l) and (1451 ), respectively. It is worth noting that Z\ ■ ■ ■ z^ and /x* satisfying both (fl"8l) and (IT9T ) are 
uniquely determined since 7V+ 1 variables are solutions of N + 1 independent equations and In (1 + -Zi) — 
Y^- is a monotonically increasing function of Zi. In addition, since 1 — Tq = X] r «* from (1461 ) and r* = Tq 



*7i 



from (1471 ). it follows that 

from which we obtain ([Tot . Finally, we obtain (fTTT ) from (|47T ) and (|48l) . 
This thus completes the proof of Proposition 14.11 
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